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DATA  REDUCTION  FOR  THE  FREE  FLIGHT  SPARK  RANGES 


ABSTRACT 

The  data  reduction  process  for  the  Free  Flight  Spark  Ranges  is  described 
with  emphasis  on  recent  modifications.  The  most  important  modification  which 
is  that  of  swerve  reduction  is  treated  in  some  detail.  Criteria  for  quality 
of  results  are  discussed  and  a table  of  time  required  for  present  machine  data 
reduction  routines  is  included . 


INTRODUCTION 


One  of  the  basic  problems  .of  the  free  flight  range  technique  is  that  of 
data  reduction.  This  is  partly  due  to  the  fact  that  in  range  work  the  actual 
motion  of  a missile  in  flight  is  observed,  and  the  aerodynamic  forces  and 
moments  are  inferred  from  this  motion.  Due  to  its  importance  for  the  inter- 
pretation of  range  firings,  the  cumbersome  nature  of  the  analysis  actually 
limits  the  output  of  the  ranges.  For  this  reason  a great  effort  has  been 
made  to  refine  the  whole  process  and  to  make  as  much  use  of  large  scale  com- 
puting machines  as  possible. 

One  book  [l]  and  two  reports  [2]  , [3]  have  described  most  of  the  phases 
of  such  an  analysis  which  now  is  somewhat  out  of  date.  This  report  will 
attempt  to  describe  completely  the  present  data  reduction  with  an  emphasis  on 
points  which  were  not  covered  in  previous  publications. 

The  report  will  be  divided  into  five  sections: 

1.  Determination  of  atmospheric  conditions  and  the  geometrical  motion. 

2 . Drag  and  roll  reductions  . 

3.  Yaw  reductions. 

4.  -Swerve  reductions. 

5.  Criteria  for  quality  of  results. 


DETERMINATION  OF  ATMOSPHERIC  CONDITIONS  AND  GEOMETRICAL  MOTION 

Before  discussing  the  necessary  measurements  and  computations  we  will 
first  describe  the  two  Free  Flight  Spark  Ranges,  A more  complete  description 
may  be  found  in  [ll]  . 

The  first  range  to  be  put  into  operation  was  the  Aerodynamics  Range . 
Through  it  models  with  body  diameters  or  fin  spans  up  to  57mm  are  launched 
from  a variety  of  guns  and  are  observed  by  as  many  as  46  spark  photographic 
stations  located  over  a 285*  portion  of  the  trajectory.  Figure  1 shows  the 
range  as  seen  from  the  gun  position.  The  cylinders  on  the  left  are  spark 
boxes  which  provide  a spark  of  less  than  one  microsecond’s  duration.  The 
missiles  will  normally  pass  inside  the  brass  frames  on  the  righo  and  trigger 
the  spark  either  electrostatically  or  electromagnetically.  Photographic 
plates  placed  on  the  steel  supports  of  each  station  then  record  the  shadow 
Of  the  missile  directly,  and  indirectly  by  means  of  the  mirrors . Figure  2 
illustrates  this  arrangement.  At  up  to  twelve  stations,  the  time  intervals 
between  spark  discharge  can  be  recorded  to  a least  count  of  5/8  microsecond. 
The  stations  are  surveyed  to  a positional  accuracy  of  .01  inch  and  an  angular 
accuracy  of  3 minutes. 

The  Transonic  Range  is  shown  in  Figure  3.  It  observes  680  feet  of  tra- 
jectory, contains  25  spark  stations  and  can  launch  missiles  up  to  eight  inches 
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in  diameter.  Its  instrumentation  is  complicated  by  the  introduction  of  a 
second  spark  and  also  cameras  which  record  the  shadows  cast  on  12'  and  1$' 
beaded  motion-picture  type  screens . In  Figure  U the  geometrical  details 
of  a spark  station  are  shown.  Due  to  the  greater  distances  involved  and 
increased  light  required  the  sparks  are  about  two  microseconds  in  duration. 
Spark  interval  timing  pulses  can  be  recorded  at  sixteen  stations.  The 
surveyed  distances  in  the  Transonic  Range  have  an  accuracy  of  about  ,01' . 

An  angular  accuracy  of  6 minutes  of  arc  is  obtained, 

a.  Atmospheric  Conditions 

For  theoretical  considerations  it  is  necessary  to  know  the  velocity 
of  sound  and  density 'of  air.  It  is  attempted  to  keep  them  constant  through- 
out the  observed  trajectory  by  air  conditioning  in  the  Aerodynamics  Range  and 
by  complete  insulation  and  heating  in  the  Transonic  Range , Before  each  firing 
the  temperature,  pressure,  and.  relative  humidity  are  measured.  From  the 
temperature  it  is  then  possible  to  calculate  the  velocity  of  sound  in  dry  air, 
VSD>  by  the  relation 

VSD  3 65,75  ^ (1) 

where  Vg^  is  in  feet/sec;  T is  the  temperature  in  degrees  Kelvin. 

In  [1U]  the  following  small  correction  for  humidity  is  derived: 

Vs  - (1  ♦ .1U  -?)  VSD  (2) 

1 

where  Pw  is  pressure  of  water  vapor 
P is  total  air  pressure 

Due  to  the  presence  of  water  vapor  the  calculation  of  density  is  of  some 
interest.  If  P is  the  total  pressure,  Pa  the  partial  pressure  of  water  vapor, 
by  use  of  Dalton's  Law  of  Partial  Pressure  we  can  write: 

p * r + p 

a w 

But  if  we  assume  the  perfect  gas  law,  this  can  be  written  as 
p " pa  + pw 

where  the  densities,  p. } are  for  the  same  temperature,  T,  but  different  pres- 
sures, namely  P,  Pa,  F^t  respectively.  We  now  refer  p • , to  standard  con- 

ditions  TQ  of  0°C  and  PQ  of  atmospheric  pressure; 


1 This  pressure  is  easily  calculated  by  means  of  the  relative  humidity  and 
tables  of  water  vapor  pressure  for  saturated  air. 
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(3) 
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Equation  (3)  together  with  the  known  values  of 


pao  Tc 


and 


rwo 


ao 


= .3783 


now  allows  us  to  compute  air  density  from  measurements  of  P,  T,  and  P^. 


b.  Plate  Measurement  and  Geometrical  Calculations 

On  both  ranges  the  z axis  points  down  range  from  the  gun,  the  x axis 
to  the  left  looking  downrange  and  yaxis  up.  Each  station  has  a local  origin 
located  on  the  intersection  line  of  the  planes  of  either  the  photographic 
plates  for  Aerodynamic  Range  or  the  screens  for  the  Transonic  Range  so  that 
its  xy  plane  contains  the  sparks^-. 

11"  x lU”  plates  are  employed  on  the  Aerodynamics  Range  and  are  measured 
on  ruled  grids  set  in  light  boxes.  Figure  5 shows  part  of  a plate  including 
a fiducial  bar  with  three  nicks.  The  edge  of  the  bar  is  parallel  to  the 
z axis  and  the  x,y  plane  is  located  by  the  nicks.  Knowing  the  distance  from 
the  edge  of  the  bar  to  the  local  z axis  we  can  therefore  make  measurements  in 
the  local  coordinate  system.  The  actual  measurements  taken  are  the  location 
of  a reference  point  on  the  shadow's  axis  and  the  slope  of  the  axis  with 
respect  to  the  fiducial  bar. 

On  the  Transonic  Range  the  usual  plate  size  is  U"  x 9' . These  are 
usually  measured  pn  Mann  Comparators  (Figure  6)  or  a Telecomputing  Telereader 
(Figure'  7)  . The  Telecomputing  equipment,  which  is  operated  by  the  Measurement 
Analysis  Branch  of  the  Computing  Laboratory,  has  the  desirable  characteristic 
of  IBM  card  output.  Since  only  positions- may  be  measured  by  the  Telereader, 
the  slope  of  the  axis  is  obtained  by  the  measurement  of  two  points  on  the 
shadow's  axi3.  Figure  8 is  a sample  plate  from  the  Transonic  Ranged.  Measure- 
ments are  made  with  respect  to  the  crossed  survey  wires  and  are  converted  to 
distances  on  the  screen  by  means  of  a magnification  factor.  These  results  are 
then  transformed  to  the  coordinate  system  formed  by  the  intersection  of  ex- 
tensions of  the  screen.  For  both  ranges, therefore,  the  measurements  can  be 


1 On  the  Aerodynamics  Range  both  the  actual  spark  and  its  virtual  image  in 
the  mirror  are  considered. 

2 Since  the  camera  is  focused  on  the  screen,  the  clear  image  is  of  the  shadow 
while  the  blurred  image  is  of  the  missile  itself. 
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reduced  to  locations  of  points  on  the  shadows  of  the  missile's  axis  and  the 
slopes  these  shadows  have  with  reference  to  a station  coordinate  system. 

The  problem  is,  then,  to  derive  from  these  data  the  missile's  location  and 
orientation  in  space  at  each  spark  station. 

Since  the  sparks  are  located  in  each  station's  xy  plane,  they  have 
coordinates  (cp  Cp  0)  for  the  spark  opposite  the  vertical  plate  and 

(c^ 1 1 , 0)  for  the  other.  (On  the  Aerodynamics  Range  the  second  spark 

is  actually  a virtual  spark  which  is  located  behind  the  mirror) . If  we 
denote  the  coordinates  of  the  shadow  of  the  reference  point  (x^,  y^,  z^)  by 

(0,  yv,  Zy)  on  the  vertical  plate  and  by  (x^,  0,  z^)  on  the  horizontal  plate, 

points  on  the  line  between  the  sparks  and  respective  shadows  must  satisfy  the 
following  equations: 


x . o y “Py  a - av 

— <—»  m m - ■ ■!— w. 

cl~^  c2  " ® = 

x mxa  y - 0 a aH 

c^*-  X^  Cg‘-  0 0 ” ZJ| 


CIO 


(5) 


Since  the  reference  point  mu3t  lie  on  both  lines  we  have  the  following 
four  equations  for  the  three  coordinates  of  this  point1 


where 


* *H  * yR 

(6) 

yE  " yv  * Ry  3% 

(7) 

■a  * *h  " c£.  yR 

(8) 

■v  .. 

^ - 2 V "cj  *R 

(9) 

p Cl'  1 

“x  C2.  “ cj' 

p °2  1 

^ ■ °r  ^ 7t 

1 The  system  is  overdetermined  because  the  lines  are  restricted  by  the 
assumption  that  they  intersect  - The  survey - and  measurement  accuracy 

iriaj  be  checked  by  a comparison.;  •'  ’the  two  values  vofz^. 


8 


If  we  denote  the  direction  cosines  of  the  missile's  axis  by  (m,  n,  p) 
then  the  equations  of  planes  containing  the  missile's  axis  and  the  shadow 
of  the  reference  point  on  the  vertical  and  horizontal  planes  respectively 
are-*-: 


x " C1 

y - c2 

0 

1 

0 - 

7V  “ °2 

<r 

i 

o 

« 0 

(10) 

m 

n 

p 1 

x - C-j/ 

y - c2* 

0 

1 

H 

0 

1 

4* 

0 - Cg1 

*H  -° 

■ 0 

(11) 

m 

n 

P 

If  tan  V is  the  slope  of  the  shadow  of  the  axis  on  the  vertical  plane, 
the  point  (0,  yy  + tan  V,  Zy  + 1)  is  a point  which  satisfies  equation  (10). 

If  tan  H is  defined  similarly,  the  point  (xy  + tan  H,  0,  Zg  + 1)  lies  in  plane 

described  by  (11) • After  a few  algebraic  manipulations,  these  equations 
reduce  to 


m/p  - tan  H - n/pj  ~T  " 

fc2 

n/p  « tan  V + m/d  — - Yy/^  * 


+ tan  H 


~ tan  V 


(12) 


(13) 


2 2 2 
and  m + n + p « 1 


(1U) 


1 These  determinants  are  the  scalar  triple  products  of  the  vectors  (x  - c^, 
y - Cg,  z - 0),  (0  - Cp  yy  - c2,  Zy  - 0),  (m,  n,  p)  and  the  vectors 
(x  - cx»,  y - c2»,  z - 0),  (xg  - c1»,  0 - Cg ' , Zg  - 0),  (m,  n,  p)  re- 
spectively. Since  scalar  triple  products  represent  the  volumes  of  the 
boxes  formed  with  the  three  vectors,  these  quantities  will  be  zero  only 
when  the  point  (x,  y,  z)  lies  in  the  plane  determined  by  the  vector 
(m,  n,  p)  and  either  the  vector  (0  - c^,  yv  - c2,  Zy  - 0)  or  the  vector 

(xg  - c^',  0 - c2‘,  Zg  - 0). 
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From  (12)  - (lU)  m,  n,  and  p may  be  determined.  If  L denotes  the  distance 
between  the  reference  point  and  the  center  of  mass  (xcm,  ycm,  zcrn),  then 

xcm  ’ % * mL  * xi 

" yR  * ^ * n (15> 

zcm  * ZR  * PL  + 2i 

where  (x..,  y^,  z.)  are  the  coordinates  of  the  local  origin  of  the  ith  station 

relative  to  a fixed  origin.  The  computations  described  by  equations  (6)  - (l£) 
have  been  coded  for  the  Crdvac  and  the  IBM  CFG  for  measurements  obtained  by 
Viann  Comparator  or  light  box . At  the  present  time  the  card  output  of  the 
Telereader  can  be  processed  by  only  the  Ordvac . 

Finally  It  is  necessary  to  compute  the  two  components  of  the  yaw  from  the 
direction  cosines . The  yaw  angle  measured  in  range  work  is  defined  to  be  the 
angle  from  the  tangent  to  the  trajectory  to  the  missile’s  axis.  For  the  flat 
trajectories  encountered  in  spark  range  work  and,  for  small  yaws  this  can  be  done 
by  the  use  of  the  simple  relations: 

XH  * a - i (16) 

\ - » - H <i7> 

where  X^  is  horizontal  component  of  yaw  in  radians 
Xv  is  vertical  component  of  yaw  in  radians 

is  "horizontal  slope"  of  trajectory 

is  "vertical  sloj  of  trajectory 

For  hand  reductions  ^ and  ^ may  be  computed  by  differencing  values  of  x and 

y between  successive  stations.  As  coded  .for  the  Ordvac  the  stations  are  divided  into 

four  or  five  groups,  and  for  each  group  x and  y are  fitted  by  a quadratic 

function  of  z.  .From  these  resulting  equations  -j£  and  may  be  obtained. 
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DRAG  AND  ROLL  REDUCTIONS1 


a.  Drag  Reduction 

2 2 

The  drag  force  as  measured  along  the  trajectory  is  defined  to  be  pd  u^Kjj, 
where  p is  the  air  density,  a the  missile's  diameter,  u^  the  component  of  the 
missile's  velocity,  unresolved  along  its  axis  of  symmetry,  and  its  drag 
coefficient.  If  we  assume  a flat  trajectory  directed  along  the  z axis, 

Ui  i u - |§  for  small  yawing  motion,  and 


. 2 2 
mu  =-pd  uK_ 

If  we  define  u'  as  ^ tfhere  p = ^ is  non-  dimensional 


and 


3p  *nere  p = 3 

“-tin 


(18) 


(19) 

(20) 


Since  distance,  not  time,  is  a more  fundamental  and  convenient  measure- 
ment in  range  work,  the  distance  is  taken  to  be  the  independent  variable  * It 
has  been  found  that  a cubic  equation  in  distance  fits  the  timing  data  quite 
well  and  for  this  equation  we  have 


t ° aQ  + ax  p + agp2  + a^p3 
dz 

„ dz  _ 3p  d 

u ’ 5E  W *V 
3p 


+ 2 agp  + 3 


pdJ 

( A) 


2a0  + 


6a3  p 


pd-’  a1  + 3a-p 


(21) 


(22) 


(23) 


1 The  drag 
detailed 


2 K's  will 


reduction  is  also  described  in  [l]  , [2]  , [3] 
description  of  the  roll  reduction  is  in  [10]  . 


always  be  related  to  J 


i 

s 


by  Ki 
m x 


J.  . 

x 


and  [U]  while  a 


11 


From  equations  (22),  (23)  and  the  velocity  of  sound  we  can  compute  the 
o Mach  number  and  its  associated  drag  coefficient.  For  the  greatest  accuracy 
we  normally  evaluate  both  at  the  center  of  the  observed  trajectory. 

Equation  (21)  is  written  on  the  assumption  that  t can  be  represented 
by  a cubic  in  z.  It  is  known,  however,  that  is  a function  of  6^,  the 

squared  magnitude  of  the  yaw,  and  M,  the  Mach  number.  Although  the  Mach 
number  variation  on  the  range  is  usually  small,  the  yaw  can  undergo  large 
changes.  Therefore,  a yaw  drag  reduction  was  formulated  by  E.  J.  McShane  [lj. 
We  will  describe  a modification  of  this  reduction.  First  the  assumption  is 
made  that  is  linearily  dependent  on  and  M: 


\ * *2V  * 


) J. 


dm 


(2U> 


where  J_  , Jn  - and  Jn  are  constants. 

V DM 


Substituting  (2U)  in  (20)  and  integrating. 


u I Jr  dp 

-v*v  /*-U 

S 


(u  - uQ)  dp 


♦i 


(25) 


1-22  " 

We  now  replace  the  quadratic  term  by  j p where  JQ  is  an  average  value 

of  the  drag  coefficient.  Using  the  fact  that  i ■ 35  = (j),  integrating 

and  rearranging, 

-2  .»  —/>■>•  3 


‘ ■ *0 X>  h (p)  * ** (p) 


(26) 


1 


I,(p> 


r r 


o dp  dp 


^(P) 


f jf 


u - u. 


l- 


*)  dp  dp 


o o 

Assuming  that  u can  be  approximated  ty  (22), 
u - u _ . 

■ "T~  * ^ p * S 7 

*1  . *1 
2a„ 


^2p 
a2 


^(p)  j ^5 


-1) 


For  convenience  we  select  Jn  to  be  the  value  supplied  by  the  standard 
2^ 

drag  reduction,  namely  — , Therefore 

al 

OOO  0 H 00 


where  M 


*S*1 


is  the  average  Mach  number  sinoe  p is  selected,  to  be 


(26*) 


aero  at  the  middle  of  observed  trajectory. 

It  now  remains  to  evaluate  I^(-p) . 

* 

In  [ 1 ] or  [ $ ] , it  is  shown  that  for  a symmetrical  missile  the  total 
yaw  X can  be  written  as 

X = \H  + iXv  - Kx  e (<tl  + 1 * p + Kg  o ^"*2  + ^ 2 )p  ♦ XR 

where  are  complex  constants  „ 

are  real, constants 
<j[£  are  real  linear  functions  of  p 

X^  is  the  "yaw  of  repose"  and  is  determined  by  a yaw  reduction  as 
described  in  the  next  section 

If  Xp  is  small, 

52  - X X - K1q2  e~2alp  + K2q2  e “2a2P  + ^(“l  + «2>P 

cos  [(gr2'^')p  + ((*20-?10^ 

whore  ■ K^Q  e ^ io 
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(2?) 


(28) 


I 


If  (28)  is  placed  in  the  definition  of  1^  and  the  indicated  integrations 
are  performed,  there  results^-: 


I^p)  * 


5-2giP  -1  + 

{2^f 


2alP  + 


r2a2p 


(2cu)‘ 


where  r 


+ + °2^P  cos  [(^2*  " $i')p  + yJ  “ cos  y 

- p r cos  (y  +^\  ) [ ( 

" ^l')2  + K + °2)2 


= arc  cos 


-(«!  + °2A 


- arc  sin 


” 2 1 


The  procedure  for  yaw  drag  reduction  is_first  to  perform  a standard 
drag  reduction  in  order  to  obtain  values  of  jjj  and  M which  together  with 
an  estimated  Jn  determine  the  cubic  coefficient.  A least  squares  is  then 

M d d d 
min  for  the  remaining  unknown  coefficients,  tQ,  — , — - J_  , and  — , 

uo  uo  uo  u62 

If  J is  not  fixed,  a least  squares  fit  of  five  parameters  can  be  attempted. 

M 2 
This  latter  procedure  has  been  ineffective  up  to  the  present  time  . 

The  variation  of  drag  with  yawing  motion  for  a group  of  identically 
shaped  models  fired  at  about  the  same  Mach  number  can  be  determined  rather 
simply,  however,,  from  the  standard  drag  reduction.  We  assume  that  their 
individual  drag  coefficients  are  average  J^’s,  over  the  length  of  observed 
trajectory.  Eq.  (2U)  can  be  written  in  the' rorm 


J-  ■ J_  +6  J n 
D D.  Dc2 

O o 


(2U‘) 


1 The  cosine  term  is  integrated  as  though  0'.  were  constants. 


2 This  yaw  drag  reduction  has  been  coded  for  the  Bell  Relay  Computer  by 
J.  Schmidt  and  L.  Schmidt. 


U 


where  8 


Pl  -Pf 


q dp  (mean  squared  yaw) 


PF 


Pp  is  p coordinate  of  missile  at  the  first  timing  station 
is  p coordinate  of  missile  at  the  last  timing  station 

— * 

The  value  of  6^  is  a by-product  of  the  first  integration  involved  in 
obtaining  Eq.  (29): 


7 


1 


2 /e20!11?  - e2glPL'  ^ „ 2 ,e2o2pF  - i2a2pL 


KW  < 


) ♦ P 


0 


where  the  integral  of  the  cosine  term  is  neglected  for  the  usual  large 
values  of  p^  - p-p  encountered . The  drag  coefficients  for  each  round  are 

then  plotted  against  the  corresponding  values  of  the  mean  squared  yaw.  A 
line  is  fitted  to  these  points  and  good  determination  of  and  K~  _ for 
the  common  Mach  numbers  can  be  made.  o 

b.  Roll  Reduction 

The  equation  of  rolling  motion  for  a missile  with  rotational  symmetry 
may  be  written  in  the  form  [10]  : 

v*  - Dv  + k^2  JA  « D v 

co^d 

"u‘" 


where  v » 0 ' 


(30) 


9 is  roll  angle 

axial  angular  velocity  _2 

TA 


-2  A ^1  ^An 

JD  - kj.  JAJ  D - D ♦ 


-2  _ md 


-jr-  (k^  is  axial  radius  of  gyration  in  calibers) 


m = mass 

A * axial  moment  of  inertia 


and  Ka  and  KA  are  defined  by  the  relation:  axial  aerodynamic  moment 

*3„2  rK  . vK  ■ 

Ao  A 


pdV 
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(31) 


The  solution  of  ( 30)  can  be  easily  obtained  and  is 
0 - B + sp  + Ae1^ 

kl  JA 

where  s ^ — 0 (steady  state  roll  per  caliber)  . 

For  range  work  the  roll  angle  may  be  measured  by  means  of  two  pins  of  different 
shape,  in  the  base  of  the  missile  (see  Fig.  8).  The  location  of  these  pins  in 
space  may  be  calculated  by  use  of  a slightly  modified  center  of  mass  reduction 
coding  and  the  orientation  of  the  vector  between  the  pins,  then  provides  the 
roll  angle. 

If  D is  sufficiently  small'*'  the  exponential  may  be  expanded  as  a cubic  in 
p and  the  standard  drag  reduction  coding  which  is  available  both  on  the  Ordvac 
and  the  Bell  Computer  can  be  used  to  fit  9 as  a function  of  z.  In  general, 
however,  a new  reduction  procedure  is  required.  Furthermore  since  the  unknown 
coefficient  D appears  in  a non-linear  fashion  the  iterative  method  of  differ- 
ential corrections  is  needed. 

In  order  to  start  the  iterations  a set  of  initial  values  is  a prerequisite »' 
To  do  this  we  differentiate  (31)  and  use  the  result  to  eliminate  the  exponential: 

- 9 * (B  - jj)  + sp  + p 0 (32) 

i 

9 can  be  computed  by  numerical  differentiation  using  first  differences  and 
equation  (32)  fitted  by  the  routine  least  squares  since  it  is  linear  in  un- 
knowns (B  - jj),  s,  and  jj.  From  these  B,  s,  D follow  and  by  use  of  (31)  for 

a particular  station,  A is  determined. 

From  equation  (31)  a relation  in  the  differential  corrections  of  these 
values  can  be  derived: 

A9  ■ 9 observed  * 9 computed  - A B * pAs  4 eDP  A A * Ap  eDP  AD  (33) 

9 is  computed  from  (31)  by  use  of  the  initial  computed  values  of  A,  B,  s,  D or 
the  values  obtained  from  the  preceding  iteration.  Since  (33)  is  linear  in  the 
correction,  ordinary  least  squares  apply  and  a set  of  corrected  coefficients  is 
obtained.  If  the  initial  values  are  close  enough,  the  process  will  converge. 

Most  rounds  require  no  more  than  two  iterations . The  c omplete  roll  reduction 
is  coded  for  the  Ordvac  while  the  iterative  process  of  (33)  is  coded  for  the 
Bell  Computer. 


1 This  is  the  case  for  bodies  of  revolution. 
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In  [l6]  the  variation  of  and  with  Mach  number  considered. 

o 

There  it  is  shown  that  can  be  quite  well  approximated  by  an  inverse 

linear  function  of  Mach  number.  For  simplicity  we  will  make  the  quite 
reasonable  approximation  that  B itself  is  such  a function. 

. . D » (a  + bM)  ^ 

V- 

Nor  from  Sq,  (22)  it  is  easy  to  show  that  a good  approximate  relation  be- 
tween Mach  number,  M,  and  position  on  range,  p,  is 

M = Mo(l  - JDp) 

Fran  these  equations  we  have 


1 + YP 

where  DQ  * (a  + bM0)“1  (value  of  D at  p ■ 0) 


-Jr 


1 + 


bM, 


-2 

The  further  assumption  that  the  ratio  of  k,  J.  to  D is  a linear  function 

i aq 

of  M or  p must  now  be  made . 


-2 


\ JA  - - D(sq  + s'p) 


Fran  these  assumptions  the  following  revised  Eq.  (30)  results 

’o  - s'p) 

( V-  s - s'p) 


v ' ■ D(  v - sQ  - s'p) 


1 - yP 

This  can  be  integrated  to 

9 - B + sp  + Cp2  + Aef^ 

where  f (p)  ■ ^*^0  ln(l  + Y p)  * Ep(l  - XjE  * ££ . .) 


(301 


(3l') 
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I 


I 


♦i 


S * S + m 

o D - v 

r\  * 


M-A) 


Y + 


By  use  of  Eq.  (31’)  we  can  now  modify  Eqs.  (32)  and  (33). 

9 « (B  - jy)  + (s  - jy)p  + Cp2  + jy  9 * 

where  f » ■ E |~1  - yp  + (y  p)^  - ( Y p)^  • .•"]  ■ E 

*-  J 2 2 
A9  =AB  + pAs  + p2A  C + ef^AA  +A.ef^PV(l  - r».) AE 


(32 1 ) 


-Aef(p)Ep2[|  - 


AiT 


(33') 


< I 


The  procedure  is  quite  similar  to  the  regular  roll  reduction.  First  B, 
s,  C , and  E are  calculated  from  Eq.  (32*)  by  a least  squares  fit.  These  are 
placed  in  Eq.  (31*)  together  with  0 and  p at  the  station  for  which  p ® 0 
and  B is  computed.  Next  y is  either  estimated  or  computed  from  Eq.  (31')  for 
a station  at  which  p is  large.  Finally  Eq.  (33')  is  ■employed  in  as  many 
iterations  as  necessary  for  complete  convergence.  The  problem  as  it  is  being 
coded  for  the  ORDVAC  will  have  the  option  of  fixing  y (ay  « 0)  or  allowing 
it  to  vary.  This  reduction  then  will  b©  either  a five  or  six  unknowns  problem. 
Recent  experience  with  small  light  finned  missiles  definitely  shows  the  need  for 
this  modified  reduction. 


YAW  REDUCTIONS1 

For  missiles  possessing  angles  of  rotational  symmetry  less  than  1’80°  and 
planes  of  mirror  symmetry,  the  definition  of  the  linearized  force  and  moment 
system  is  well  known  [l]  , [9]  . If  the  requirement  of  mirror  symmetry  is 
relaxed  slightly^  we  have  the  spin  producing  moment  J,  Oof  equation  (30),  [10] 
If  in  addition  to  this  we  assume  either  a slight  configurational  asymmetry  or 
slight  mass  asymmetry  [13]  , the  transverse  force  and  moment  definitions 
receive  a constant  missile-attached  force  and  moment  increment.  In  the  theo- 
retical development  of  [9]]  the  usual  coordinate  system'  is  orientated  so  that 
the  one  axis  lies  along  the  missile's  axis  of  symmetry,  the  two  axis  is  in 
the  horizontal  plane  pointing  to  the  right  and  the  three  axis  is  determined 
by  the  right  hand  rule.  This  selection  of  the  2 and  3 axes  is  the  exact  re- 
verse of  the  orientation  of  the  ranges*  x and  y axes.  Since,  however,  the  yaw 


In  order  to  understand  fully  the  next  two  sections,  they  should  be  read  in 
conjunction  with  [9]  • 

2 

In  other  words  the  assumption  is  made  that  the  differential  canting  of  . 
fins  of  a finned  missile  affects  only  the  axial  moment* 
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as  defined  in  [”9]  is  measured  from  the  missile's  axis  to  the  tangent  to 
the  trajectory,  the  reverse  of  the  range  definition',  we  have  X9  = X and 

X^  « Xy.  where  Xg,  X^  are  the  yaw  components  of  |_  9 J #.  The  aerodynamics 

force  and  moment  acting  perpendicular  to  the  missile's  axis  can  now  be 
defined  as: 


(-Kjj  + ivKp)X  + (vK^p  + iKg)^  + KN  Xe  e 


id 


where  9 


F2  + 1P3  = pd2ux2 

M2  + i&3  - Pd3^2  ^(-vKj  - iK^X  + (-1^  + ivK^p  + i^  Xg  e 
,P 


(3W 


30 


r 

« 


vdp  + 9 (roll  angle  with  respect  to  force  due  to  asymmetry) 


v ■ 


*1 


(cop  w^)  angular  vel'ocity  vector 
X * Xg  + iX^  (complex  yaw) 
co.d  w_d 

|jl  =>  + i (non-dimensional  angular  velocity) 


1 


"n 

,KF 

h 

hF 

*T 

h 

Kjjt 


magnitude  of  asymmetry  angle 

Normal  force  coefficient 

Magnus  force  coefficient 

Damping  force  coefficient 

Magnus  force  coefficient  due  to  cross  spin 

Magnus  moment  coefficient 

Overturning  (or  righting)  moment  coefficient 

Damping  moment  coefficient 

Magnus  moment  coefficient  due  to  cross  spin 

Normal  force  coefficient  due  to  small  asymmetry 


Kjj  Moment  coefficient  due  to  small  asymmetry3' 


In  the  case  of  mass  asymmetry  only  the  product  X£  has  meaning  and 

can  be  computed  from  the  physical  measurements  of  the  projectile.  This 
remark,  of  course,  also  applies  to  the  product  Kf7~Xg  * 
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0 


It  Is  easy  to  see  that  the  exponential  coefficients  of  and  perform 

the  task  of  insuring  that  the  force  and  moment  associated  with  the  asymmetry 
be  missile-attached. 

The  equations  of  motion  which  can  be  obtained  from  the  above  definitions 
for  flat  trajectories  are: 


x"  + (H  - i vj  X*  + [ - M - i v t]\  - G + J€  e^ 


v » D ,v 


where 


H = JL  “ JD  + k2  ‘ JH  1 JL  " JN  " JD 


B transverse  moment  of  inertia 
M - k‘2[jM  * v2  k,2  Jp]  i V* 


T « JT  - 


V JT 


^ [JB 


+ k2  JE~±V. 


-2  md2 

kg-  * -g-  (kg  is  transverse  radius  of  gyration  in  calibers)# 


- The  solutions  for  these  equations  of  motion  can  be  written  for  the  case 
of  slowly  changing  spin  as  ? 


X - e’V  + i^10  p + ^1  jp  ^ + Xg  e "a2p  + i^20  p + ^2  2 ) +XR+K3e19  (36) 


where 


c #>■«[{-*  -< > ■ • [-  • * <**‘3  i» 

_ P r •> 

-OgP  + K^qP  + $2  1/2  J • -H  + iv  - D «1  - [-in  + 2iv  (2T-H-S)]  1/2jdp 


1 See  Eqs.  (8)  and  (11)  of  [9]  • For  the  relation  defining  K~  see  [13]  • 


J 


*R  * [«  - ' V®  JW  * * H> 


■ m ^ Jj*  1 2(2T-H-P)2] 
1 m2  + Uv2  (2T-H-B)2 


in  - v2  - UM  - H2  . 

and  Kg  are  complex  constants  depending  on  initial  conditions 

K,  * (= ; ■-  ~r ; . - K,n  e i^30"9o^;  K,„  real. 

3 * it^  -v)]  [-Oj  * i(s^2  -v)]  30  30 

If  we  differentiate  and  add  the  definitions  for  the  exponents  in  Eq.  (36) 
and  equate  the  real  and  imaginary  parts, 1 

v - ^ + $2*  where  * Ko  + ^i”  p (37) 

II  « (0l  + og)  - D e1  (38) 

Ilext  we  differentiate  and  subtract  the  definitions 
[~m  + 2iv  (2I-H-D)]  1/2  = -(c^  - Og)  + i(^  - ^) 

, ^ in  ■ ({^  ■”  ^g)2  " (a-j  - ag)2  (39) 

.\  v (2T-H-D)  - -(^  - ctg)  (^  - 0g)  (UO) 

Inserting  the  definition  of  m in  (39)  and  using  (37)  and  (38), 

.1  ,t  A ^1  A.  ... 

M ■ * (?2  “ + D(a^+Og«  -tjt  D)»  (Ul) 

Equation  (UO)  can  be  solved  for  T: 


T - -1/2  - dj)  % % - H-D 

h *4  J 


Although,  as  we  shall  see,  D can  be  computed  from  the  yaw  reductions,  the 
resulting  values  an;  sometimes  quite  bad  due  to  experimental  limitations . 
For  use  in  Eqs,  (38),  (Ul)  and  (U2),  8 should  either  be  estimated  or 
directly  measured  by  a roll  reduction. 
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It  remains  to  compute  £ ^ from  its  definition.  But  this  can  be  done 
by  means  of  equations  (36),  (39)  and  (U.0) . 


<sq'  * it,')2 

1 1 

A 

- < >2  - - V2. 

+ 2(0^  - - $2  ) 

*1  " 

r 1 i 1 

• f 

2 

1 

<r* 

+ 

CM 

/ V 

1 

1 r_i 

2 

$1  + 

lA^'J 


(U3) 


$2  are  quantities  which  are  obtained  in  the  curve  fitting 


‘process  to  be  described  shortly- » By  means  of  equations  (3?).  (38),  (Ul), 

(h2)j  and  (U3) , the  spin  v,  and  three  combinations  of  aerodynamic'  coefficients 
M,  H,  and  T may  be  calculated.  Although  can  be  directly  obtained  from  M, 

an  independent  'evaluation  of  is  required  in  order  to  obtain  and  from 

H and  T respectively.  The  variation  of  the  turning  rates  is  caused  physically 
by  the  change  in  v and  M during  flight.  Using  equations  (30),  (37),  and  (Ul) 
we  have 

a ^ + ^ 

d - I - (Wt) 

11  + 02 


M' 


. » w »«  .n 

$-1  • + & 


”2  T ^2  * ^1  ’ 
Symmetrical  (Epicyclic)  Yaw  Reduction 


¥e  will  first  discuss  the  reduction  process  for  a symmetrical  missile. 
<%>  - * Ko  * 0) . This  yields  a rosette  or  epicyclic  type  yawing 

£ £ 

motion.  Fig.  9 shows  the  yawing  motion  of  a spinning  body  of  revolution 
with  negligible  yaw  of  repose  and  Fig.  10  displays  that  of  a rapidly  spinning 
finned  missile. 

The  first  step  in  this  reduction  is  the  calculation  of  the  yaw  of  repose . 
For  velocities  greater  than  600  ft/sec  A,  a good  approximation  from  its 

definition  can  be  obtained  from  equation  (36)  and  is  XR  = . All  the 

“ Mu 

quantities  in  the  above  relation  are  either  known  quantities  or  may  be  quite 
easily  estimated.  Of  course  once  initial  values  of  ^ are  obtained,  equations 

(37)  and  (Ul),  neglecting  the  effect  of  a, may  be  used.  The  yaw  of  repose  is 
real  and  negative  for  positive  spins  which  means  that  the  missiles  will  tend 
to  point  slightly  to  the  right  of  the  trajectory .2 

1 flange  firings  are  usually  restricted  to  velocities  greater  than  this  value, 

2 It  should  be  remembered  that  the  positive  real  axis  has  been  selected  to 
point  to  the  left. 


' ZZ 


This  computed  value  of  XR  is  then  subtracted  either  numerically  or 

graphically  from  the  measured  X^'s  and  the  resulting  yaw  values  plotted  in 

groups  of  five  or  six  stations.  For  these  groups, which  cover  about  1$%  of 
the  total  length  of  observed  flight,  the  effects  of  damping  and  change  in 
turning  rates  are  assumed  to  be  negligible,  ^ * 0 ) 11518  points 

should  then  satisfy  equations  of  the  form 

#10  * <*10  (P-J>o3  t „ ^20*  <*10  (p-pJ 


where  K. 


X " XR  * ^0 
i (j, 

v *io 

■ Kio  e 5 


e 


K20 


K.  real 
xo 


pQ  is  value  of  p for  middle  station  of  group.  Multiplying  by  e 


-itf'  (p-p  ) 

° (X  - XR) 


K10e 


#10  * Wio-^'XP-Pol 


*20 


(U6) 

-¥’(p-p0) 

$20  + ^20-^  (p“Pol 

a?) 


This  can  be  done  quite  easily  by  use  of  a compass  since  this  multi- 
plication is  equivalent  to  rotation  of  the  yaw  angle  at  p through  the  angle 
-#f'  (P“PQ) • Should  be  so  guessed  that  it  is  quite  close  to  the  local 

slow  rate  ^2' q above  equation  reduces  to  the  equation  of  a circle: 


-#5o  (p-po> 


(x  - XR) 


lo 


#10  * «*10  - ^0)(p-pol  „ # 


+ K, 


20 


20 


(U8) 


1 

This  trial  and  error  process  can  usually  be  done  quite  quickly  and  the 
completed  solution  is  shown  in  Fig.  11 « Note  that  5^j*o  can  o5)^ai-nec5 

from  the 'position  of  the  points  on  the  circles.  If  the  ^Q's  obtained  for 

adjacent  groups  with  due  care  for  multiples  of  360  are  considered,  some- 
what more  accurate  values  of  *s  may  be  computed  for  the  midpoint  of  the 

interval  between  them.  If  these  values  for  all  of  the  groups  are  plotted 
against  p,  estimates  of  0?  can  be  made  by  dse  of  a straight  line  fit. 

Finally  the  Kio's  are  plotted  on  semi -log  paper  against  'p  and  can  be 

determined  from  the  slope  of  a fitted  straight  line. 


From  the  i above  brief  description,  it  can  be  seen  that  first  values  of 
the  four  initial  constants  Kio,  9^0's  for  the  range  origin,  and  the  six 

parameters  a^,-  and  (jS^  may  be  determined.  It  is  now  necessary  to  apply 


1 


An  interesting  variation  of  this  process  which  makes  use  of  an  analogue 
computer  is  described  in  . 


23 


the  method  of  differential  corrections  which  was  first  described  for  the 
roll  reduction.  The  equation  for  the  differential  correction  and  the  re- 
sulting normal  equations  are  given  in  [3j  and  [8]  .1  This  ten  unknown 
yaw  reduction  has  been  coded  for  the  Bell  Relay  Computer,  Eniac,  Edvac  and 
Ordvac , 

If  a roll  reduction  has  been^done,  much  additional  information  is  avail- 
able for  the  yaw  reduction.  The  D which  appears  in  equations  (38),  (ill)  and 
(il2)  is  always  computed  from  the  drag  and  roll  reductions  or  estimated  from 
experience  with  the  configuration  under  consideration  and  is  not  obtained 
from  Eq.  (UU) • The  two  values  are  compared  in  order  to  give  an  indication 
of  the  quality  of  the  yaw  reduction.  Insertion  of  known  values  of  8 and  v 
would  eliminate  one  {2(7  and  one  (jfV  . The  requirement  that  M'  = 0 would  then 

eliminate  the  other  . Finally  an  inspection  of  equation  (36)  shows  that 

* a,  when  v - 0.  Thus  we  see  that  the  number  of  unknowns  can  vary  from 

six  to  ten.  As  a result  of, this  fact  three  variations  of  the  basic  ten 
unknowns  yaw  reduction  have  been  coded  for  the  Bell  Relay  Computer.  ■ They 
are: 

2 , ti 

1.  eight  unknowns  reduction  which  uses  fixed  values  of  0.  as  obtained 

) 

from  equations  (Ui)  and  (U5)  and  initial  values  of  ^ . 

2 . seven  unknowns  reduction  which  uses  fixed  values  of  and  the 

additional  restriction  that  (jS^  + 3 v.is  known. 

3 . six  unknowns  reduction  which  is  based  on  the  assumption  of  zero 
spin  and  constant  moment  coefficient  (M  « 0) , For  this  reduction 

^ + #2  55  °>  \m  a2>  and  » ff”  ■ - 0 . 

For  certain  rounds  these  special  reductions  have  proven  to  be  superior  to 
the  basic  ten  unknowns. 


1 [3]  uses  the  triangular  square  root  method  of  matrix  inversion  while  [8] 
advocates  two  methods  which  make  use  of  the  extensive  symmetry  of  the  normal 
equations.  For  considerations  of  machine  time  we  have  found  that  most  of 
the  time  of  computation  is  devoted  to  input,  output,  and  formation  of  normal 
equations.  As  a result  of  this  experience  most  of  our  recent  work  has  used 
essentially  the  Gaussian  elimination  method  of  inversion.  It  has  also  been 
found  that  a floating  binary  point  routine  as  described  in  [12 J is  of  great 
value  and  does  not  greatly  increase  machine  time  requirements . 

2 When  a variation  of  the  eight  unknowns  yaw  reduction  was  coded  for  the  CPC 
at  the  Naval  Ordnance  Laboratory,  use  was  made  of  the  symmetry  by  means  of 
performing  all  operations  in  the  arithmetic  of  complex  numbers.  This  was 
possible  because  this  reduction  can  be  considered  as  a problem  in  four 
complex  unknowns. 
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Asy:  metrical  (Tricycle)  i'a w Reduction 

Since  it  is  extremely  difficult  to  reduce  the  yawing  notion  of  an 
asymmetrical  missile  without  knowledge  of  the  rolling  notion,  we  will  "only 
consider  the  reduction  of  rounds  for  which  a roll  reduction  has  been  per- 
form'd. It  should  be  emphasized  that  equation  (36)  is  exactly  correct 
only  when  v is  constant.  It  is , assumed  that  for  small  variations  of  v , 
when  it  does  not  get  close  to  $2^,  the  relation  is  Still  valid.  When  v 

varies  near  resonance  occurs  ?nd  equation  (36) < clearly  is  not  valid.1 

From  the  roll  reduction  we  knew  the  change  in  ^ but  not  The  problem 

then  is  to  compute  (jf^o  and  ^30*  | ( 

I ■ 

I f"i 

First  all  of  the  yaw  points  arc  plotted  and  rotated  through  angle  -(0-8  ) 
This  yields^  • 1 


-(e-eo)i 


(-^  + i^)P  - i(0-eo)  (-a,  + i5f2)p-i(0-eo)  i{/ 

h 0 + x2e  +Kdne 


(U?) 


♦ 

Since  0 i3  linear  or  very  nearly  linear,  the  result  is  an  epicyclic  motion 

i*C  , !h 

about  the  point  e (Figures  12,  13,  lu,  illustrate  this  process).  It 

is  now  easy  to  determine  K^q  and  since  they  are  the  polar  coordinates  of 

the  center  of  this  epicyclic  motion.  The  origin  is  moved  to  this  point  and 
the  points  rotated  back  through  the  angle  0 - 0Q.  lit  is  now  possible  to  per- 
form the  usual  epicyclic  reduction  and  so  obtain  the  remaining  parameters. 

The  tricyclic  yaw  differential  corrections  routine  has  been  coded  in  two 
forms  on  the  Doll  Relay  Computer  and  is  now  being  coded  for  the  Ordvac.  The 
reduction  as  coded  for  the  Ordvac  will  be  a nine  unknowns  problem  formed  by 
combining  the  seven  unknowns  epicyclic  yaw  reduction  with  the  two  parameters 
K30,  fly 0 °£  the  tricycle  arm.  The  two  Bell  Relay  Computer  versions  are  eight 

and  ten  unknowns  problems  and  were  formed  by  combining  the  two  tricycle  para- 
meters wi+h  t<bn  six  unknowns  and  eight  unknowns  epicyclic  yaw  reductions 
respectively. 


In  order  to  reduce  recent  firings  bv  J.  D.  Wicolaide3  and  L.  C.  MacAlliotor 
for  asymmetrical  rounds  which  go  through  resonance,  use  is  being  mado  of  the 
Exterior  Ballistics  Laboratory's  Analog  Computer  in  a fashion  similar  to 
that  described  in  [?j  - Work  is  bci rig  done  on  this  by  J.  Schmidt, 

2 

This  method  was  suggested  by  I..  C.  FacAllister  [15]  * X.  for  most  tricyclic 
models  is  oxtronoly  small  and  so  is  nog'.  cebo  i . 


SWERVE  REDUCTIONS  * j 

The  last  bit  of  dynamic  information  which  can  be  obtained  from  a free 
flight  range  firing  lies  in  the  lateral  motion  of  the  missile.  The  equation 
of  motion,  which  the  x and  y coordinates  of  the  center  of  mass  of  the  missile 
must  satisfy,  may  be  written  in  the  form: 

m(X  + iy")  - (Fx  + iFy)  - img  - mac  (50) 

where  F is  x component  of  the  aerodynamic  force 

F is  y component  of  the  aerodynamic  force 
«/ 

g is  acceleration  due  to  gravity 
ac  is  the  Coriolis  acceleration 

If  we  resolve  the  aerodynamic  force  along  the  x and  yfax&3, remembering 
that  the  small  angle  between  the  1 axis  and  z axis  is  A + --■+  ^ and  that  the 
2-3  ax6s  are  the  negatives  of  the  xy  axes  when  this  angle  is  zero: 

Fx  + iFy  =*  -(F2  + iF3)  + F1  (A  + (5l) 

But  F^  is  defined  by  the  relations: 

Fi  “ mmh  ui  “ u *da  : 

i 

In  [l]  and  [9j  the  following  equation  is  derived  for  the  symmetrical 
missile: 

A'  - in  * (-JL  + ivJF)  A + (vJyj,  + iJg)  n + y * 0 (5f3)  ; 

where  y is  a small  contribution  of  gravity. 

10 

The  insertion  of  the  term  J„  A€  e in  Eq.  (53)  for  the  asymmetric  missile 

£ , 

does  not  effect  the  approximation  A - in  * 0. 

Using  relations  (3h)  (52)  and  (53)  and  converting  from  K *s  to  J »s  where 
necessary  we  can  rewrite  (5l)  in  the  form:  ^ * 

2 r n • 

FX  + ^y  " ^ (JL  ’ ivJF)X  + + ivJXF)X’“  JN£X€  e ~ (x  + iy)  (5U)  j 

i; 

j 

1 . o l* 

The  standard  value  of  g at  the  Aberdeen  Proving  Ground  is  32, 152  ft/sec . . !; 

. ! 
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The  usual  change  of  independent  variable  is  now  made 


cn  + iyti  * (x  + iy)  * | |- 


(x  + iy)  (|)2  - (x  + iy)  £ (|)2 


x + iy  - (^)  (x"  + iy")  + (x  + iy)  ^ 

Substituting  (5U)  and  (56)  in  (50), 


c11  + iy” 
d 


3.  4 

(JL  - ivJp)X  + (-Js  + ivJ^X*  - JN  Xe  e 30  - - -5- 


1 

Integrating 


x +iy  x ® +5y  ® ^ P 

*2^  - ° d ° ♦ - -d  V|  j (JL  - ivJp)X  dp  dp  t f (-Js+ivJjp)(X-X0)dp 


* \Xe 


O J 

J v 


8 dp  dp  - - 


+ ^c} 

d d 


P P 

2 f -2 

where  yQ  ■ gd  J u dp  dp 


. J J 

* * 

4. 


agd2u~2  dp  dp 


If  we  assume  that  the  drag  coefficient  is  constant, 
u"2  «=  uo"2  e2JDp 

. ,2.  -2  e2JDp  “ 2JDp  ~1  • _j2  -2  f 2 J„p3  J„2p^ 

•y°'Bdu°  — — 'ed  ° 


The  variation  of  v in  the  product  vJ.^,  was  neglected. 
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It  now  remains  to  compute  the  Coriolis  acceleration  a^ . It  is  shown 

in  most  mechanics  texts  that  the  vector  = 2 j£x  1i  where  K is  the  angular 

velocity  of  moving  axes  with  respect  to  an  inertia  set  vOf  axes  and  u is  the 
missile's  velocity  with  respect  to  the  moving  axes.  Ji: is  therefore  a vector 
directed  along  the  earth's  axis  of  rotation  with  magnitude, -A  , of  2rt  radians 
per  day  or  7.2?2  x 10 “5  rad/sec.  We  now  introduce  an  auxiliary  coordinate 
system  formed  by  the  unit  vectors  e-, , e2,  located  at  the  gun  and  so 

orientated  that  points  up  and  e^  points  north  along  the  meridian  of 

longitude . If  ^ is  the  angle  between  and  the  ranges ' z-axes  for  which 

positive  orientation  is  defined  to  be  from  ^ to"*e^  and  the  unit  vectors 

along  the  ranges'  axes  are  denoted  by  and^z  respectively  * 


e2  = ey 

"e^  = oosp  sin  & 

(61) 

„ •> , 

Now  u = u a 

z 

^1  3 ^2  x^3  * ^x  cosJ^  s^n 

(62) 

and_fl_- Aj>3  cos  9 

+ sin 

(63) 

where  9 is  the  latitude  taken  as, positive  for  the  Northern  hemisphere 
. 4 . tTg  « .2  Ax  u 

* 2 cos  ^ - *ex  sin  P ) cos  9 + ly(sin  9)j  xuez 

» 2u  Jl  ^sin  cos  9 e^.  + sin  9 (61*) 

Using  our  complex  notation* 

aC  “ 2u  A (sin  9 + i cos  9 sin  ) (65) 


■*"  For  Aberdeen  Proving  Ground  9 » 39°  26' . The  ranges  are  so  orientated 
that  is  270°  for  the  Aerodynamics  Range  and  168°  for  the  Transonic  Range  • 
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4-  iyc  - 2 fi  d2 (sirO  + i cos  6 sin  f 


P P 

rr  r> 


'0  ^0 


u"1  dp.  dp 


n u d2  (sin  9 +icos  9 sin  iff  )(p  + 

4 6 


v>>\ 

+ ir ' 


2 . V*  J*p 


(66) 


Evaluation  of  the  Turtle  term  ^ows^they  W iwolving 

Eqs . (60)  and  (66) « it  “e,.  will  first  consider  the 

S ^hfLe  "general  tricyclic  yav^  fissile. 

Epicyclic  Swerving  Motion  (K3  = K^XC  = 0) 

. /vo\ . a <? q urne  ■ 0(n  * 0,  which  xinplies 

If  we  insert  eq.  (3 6)  an  eq.  (58);  assume  ^ * 

v * * 0 and  integrate,  the  result  may  be  written  as 


, ^ + j(r  » rc  .»  raj.  . P * Qp  * hh  . <-i  + <)p 

' ♦ a*  c<-°2  * ♦ ^4^ 

where  P and  9 are  complex  constant  depending  on  initial  conditions1 


Hi  ■ Eil  * “iZ 

Ru  - an  JL  * ’“ipF  * °il  JS  * dil  JXF  ' 


2v  tf. 


io 


T72 TT2JL  + 

Wl0  *al)- 

77?  2x2  F 

Wi0  * “i  > 

°i  y 4. 

v <0 

<0  * *±  S 

W>2  2 W 

+ -a. 

^io  i 

ai2  JL  + bi2  JF  + C 

;12  JS  + di2  JXF 

<o“i 

-75 572  l 

«ic  * °i  > 

»<-“i2>  , 

(0*  * J? 

v>io  i 

the  integrals  o f Eg.  t58>  together  withX„  «r. 

absorbed  by  p and  Q* 


(67) 
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*io 

*72 T 

i'io  - “i 


JS 


V 


<<o  * 4 > 


ZF 


*d  + . 

3 


p p 

r>  r 


u 
o o 


(JL  - ivJy)  XR  dp 


P 

dp  + ^ (-Jg  + 


i-vJjy)  \R  dp  (Drift) 


If  we  insert  the  approximation  - for  XR,  neglect  the  very  snail 

Mu  J Dp 

single  integral'  and  remember  that  f.or  constant  D and  Jp,  v • vq  e and 

1 


”^T)P 

u » u0  e , the  drift  components  may  be  computed  to  be 


P P 


P P 


XD-  r 

“t^l 


sJ 
O O 


& dp  dp- -I — 

Mu"  X Hu-  u 


(D  + 2Jn)p 


dpdp.^f^Hj) 


yD.  , 
r jf 


B dp  dp 

A Mu" 


■ - J 


gdV^B 

fm?7  j 


P P 
n r 


2(D+JD)p 


gdv 


-2 


dMp.Jp— B(| 


2 D+Jn  .3 


D p" 


Mu' 


o o 


(68) 


. If  and  Jp  are  estimated,  Xp  and  y^  may  be  computed  and  subtracted 

from  x and  y respectively.  This  then  makes  Eq.  (6?)  a linear  equation  in 
eight  unknowns  to  be  fitted  to  the  measured  values  of  x and  y at  each  station. 
This  may  be  done  by  the  standard  least  squares  technique .2  Jj  and  Jj,  nay  be 

computed  from  R^,  inserted  in  the  formulas  for  Xp  and  yD  and  the  process 

iterated.  An  iteration,  however,  is  not  usually  necessary. 


-2 


1 In  almost  all  firings  where  drift  is  measurable,  D - Jp  - k2  ~ Jp . 
Use  of  this  approximation  will  simplify  (68) . 

2 In  this  calculation  the  terms  which  appear  in  the  exponentials  are 

replaced  by  + ^i  j^),  This  inconsistency  in  the  handling  of  the 

fl^'s  was  felt  to  yield  a better  approximation  to  the  true  motion.  See  page 
635  of  [1]  or  page  20  of  [2]  . 
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For  missiles  possessing  a fairly  high  rate  of  roll,  0,  is  usually  rather 
large  and  for  these  missiles  the  size  of  its  associated  swerving  motion  Ir^kJ 
is  quite  small  and  poorly  determined . But  a knowledge  of  R?  alone  puts  usJ‘ 
in  the  position  of  having  to  solve  two  equations  in  four  unknowns  1 The  situa- 
tion is  improved  by  the  fact  the  d21  and  dgg  are  quite  small  and  hence  J.^  can 

be  neglected.  Since  c21  and  c22  are  usually  small,  can  also  be  neglected. 

In  those  cases  where  c22  is  not  small  enough  to  neglect  Jg,  an  estimate  of  Jg 

is  inserted.  The  remaining  set  of  two  equations  in  two  unknowns  i3  then  solved 
to  yield  values  of  and  . 

For  missiles  which  have  been  launched  with  no  roll  an  important  simplifi- 
cation is  possible.  For  these  missiles  it  can  be  seen  from  Eqs . (3$)  that 
$1  “ ~$2  an^  °i  * °2  * additional  assumption  M 1 * 0 provides  the  result 

that  0^  * 0*  Insertion  of  these  relation^  in  the  definitions  of  R^  now 

results  in 


or 


% 


*11 


*12 


(69) 


R, 


(nf!?  _ 

vrl0 
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U'lO  + 


°1‘2)  _ 

“772  JI 


* -R, 


V) 


,»2  ; _ 2 s 


Ko + V 


(70) 
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'10 


TTTT — ? 
n.o  + °i 


By  means  of  Eq.  (69)  the  swerve  curve  fitting  of  Eq.  (67)  is  now  reduced  to 
one  involving  six  unknowns .1  Both  the  six  and  eight  unknowns  swerve  re- 
ductions have  been  coded  for  the  Bell  Relay  Computer  although  only  the  eight 
unknowns  swerve  reduction  has  been  coded  for  the  Edvac. 

Although  the  eight  unknowns  swerve  reduction  Just  described  has  been 
employed  successfully  for  some  years,  it  possesses  three  undesirable  features. 
They  are: 

1.  The  reduction  is  so-  formulated  that  the  contribution  of  the  higher 

frequency  mode  of  the  yawing  motion  to  the  swerving  motion  has  to  be  neglected. 

2.  In  order  to  correct  for  drift  it  is  necessary  to  guess  a K~  and  a K~. 

:It  would  be  desirable  to  use  drift  in  the  evaluation  of  and  Kp.  L * 

3.  The  effect  of  / 0 is  handled  in  an  inconsistent  manner. 


1 Note  also  that  for  no  spin  - y^  - 0. 
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It  was  decided  to  eliminate  features  (l)  and  (2)  by  placing  equation 
(67)  in  a form  in  which  the  aerodynamic  coefficients  themselves  and  not  the 
intermediate  E^'s  are  the  unknowns.  In  doing  this  it  was  further  decided  to 

neglect  the  affect  of  K^p . Finally  the  inconsistent  treatment  of  non-zero 

a was  corrected  by  means  of  an  asymptotic  series  expansion  of  the  swerve 

integrals . 

In  order  to  obtain  this  expansion  we  will  consider  the  integral  P e%p 
where  p"  is  constant.  If  this  integral  is  integrated  by  parts,  ^ 


J e?dP  * f <*P 

eP  f e*3 


eh' 


dp 


Pf 


P 

—3-  + 3 P"2  [4  dp  * e^V"1  + p"  p’“3)  (71) 

p*-3  ^ p*4 

e 


The  neglection  of  the  remainder  term  is  valid  for  small  *-—0  since  its  mag- 

. ...  . . p . 9.1* 


nitude  is  approximately  3 


11 

P'2 


P' 

An  examination  of  the  process  of 


generation  of  this  series  shows  that  the  magnitude  of  the  reminder ^term 
at  the  nth  step  is  approximately  1*3*5 (2n  - 1)  R ” * 

I £. 


this  we  see  that  for 


P 


i2 


J.  A 

#t"lP 


From 


< 1 the  reminder  will  first  decrease  and  then 


increase  without  bound.  This  property  is  a characteristic  of  an  asymptotic 
series . 

11 

For  the  small  0 usually  encountered,  however,  Eq.  (71)  is  a good 

P'Z 

approximation  to  the  integral.  If  we  integrate  (71)  again,  neglecting  squares 

f 


JJeVp  * 4 * ap 


0V 

P‘fc  p'^  -J  P 

■*  (p*’2  + 3p"  p’4*) 


P 


(72) 
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By  means  of  (71)  and  (72),  (£8)  can  be  written  as: 


x + xr  + i(y  + yr  + yr) 

2 a ^ - P ♦ Qp  + JL  S1  + v0  JF  S2  4 Js  s3 


where'*' 


S1  • h e1  [pi-2  + 3 P"l  Pi' 

,U]  + Kg  e^2 

[Pi"2  + 3?i  Pi4* 

1 

J 

A 

S2  - Kj  e1  [?{-2  + 3f>£  P{'U_ 

] +K2eP2 

Pi'2 

♦.a. 

V0JFd 

S3  ’ Ki  °h  [fi'1  + 'Pi  Pi'3_ 

h r 

+ K,  e c 

/r1  +pspr3] 

Pi  - -aiP  + i(^0P  + <t'l  £-  ) 

(73) 


?i  ■ (D  - a±)p  + i^ioP  + 0±  ^') 

The  definitions  of  the  S.'s  together  with  equation  (68)  show  that  they  are 
known  functions  of  p.  Eq,  (73)  is  then  a linear  equation  in  seven  unknowns 
and  can  be  fitted  by  the  usual  least  squares  method.  This  seven  unknowns 
swerve  reduction  will  be  referred  to  as  the  LFS  swerve  reduction.  For  certain 
rounds  either  Jj,  or  Jg  can  be  omitted  thereby  providing  us  with  LS  and  LF 

swerve  reductions  respectively.  The  LF  reduction  has  been  particularly  suc- 
cessful and  is  the  standard  epicyclic  swerve  reduction  coded  for  the  Ordvac. 
The  IS  swerve  which  for  zero  spin  is  equivalent  to  the  six  unknowns  swerve 
has-been  coded  for  the  Ordvac  also.  Both  the- LFS  and  LF  swerve  reduction 
have  been  programmed  for  Bell  Relay  Computer. 


Tricyclic  Swerving  Motion 

If  the  complete  Eqs . (36)  and  (58)  are  combined,  the  following  expanded 
form  of  Eq.  (67)  for  non-zero  constant  spin  may  be  obtained. 


x + XC  + *T)  * + yG  + yC  + yD^ 


- P + 


*hh 


(-®i  + i^o)p 


1 From  Eq.  (68)  we  see  that  , is  (*■*■  + — j—  pJ)  and,  hence,  is 

. La  tti0  0 

independent  of  J^.  This  term  comes  from  J'J'  dp  dp.  A similar  obser- 
vation applies  to  & . 


(-Og  + i(2f20)p 

+ hh  e + R3K3 

gdv„  J2.  D+2J. 


iB 


m 


voJFd 
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where  R^ 


JH>€  , 

-Tj-  ) ♦ i»  (Js*  Jy)  * V2 


Eq.  ( 7U)  presents  a routine  ten  unknowns  curve  fitting  problem  for  least 
squares.  This  has  been  coded  for  the  Bell  Relay  Computer.! 


For  zero  spin  Eq.  ( 7U)  is  invalid.  A second  reduction  (L-L)  has  been 
coded  for  the  Bell  Relay  Computer  for  use  on  tricyclic  rounds  with  zero  or 
near  zero  spin,  (|vj  < 10 “M  . This  reduction  was  derived  on  the  additional 
simplifying  assumption  that  all  the  aerodynamic  coefficients  other  than 

may  be  neglected.  The  near  zero  constant  spin  assumpti on. allows  the  re- 

placement  of  e in  Eq.  (36^  by 


i) 

e 


o 


ivp 

(e  ) - e 


o 
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o t 
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,b 


n 


h 1 


+ i(vp 


For  the  integration,  only  the  first  term  of  Eq.  (72)  is  used  and  the  drift 
terms  are  omitted. 

x + xu  + i(y  + yQ  + yQ) 


* P + ^ + \ S1  + jl2  Sh 


(75) 
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Since  it  was  felt  that  JT  K,  is  well  determined  by  the  swerving  motion,  both 

Eg  J 

Jj.*s  are  solved  for  in  the  hope  that  their  ratio  would  provide  a correction 


bO  * 


1 More  again  Qp  terns  appearing  in  the  exponentials  are  replaced  by 

iO/ioP  ♦ > • 


3h 


* ■ -*\**V-  ’ 
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Finally  we  will  describe  the  tricyclic  swerve  reduction  which  is  being 
coded  for  the  Ordvac . This  reduction  will  be  based  on  the  ideas  which  were 
developed  in  the  various  swerve  reductions  already  discussed.  It  has  been 
decided  to  reduce  the  four  degrees  of  freedom  expressed  by  and  Rg  of 

Eq.  (74)  to  two,  namely  and  Jg,  revise  the  tricyclic  term  so  that  zero 

spin  is  no  longer  a special  case  and  incorporate  the  full  corrections  of 
non-zero  /J 's . With  these  aims  in  mind  the  following  equation  may  be  derived, 


x + xQ  + Xjj  + i(y  + yQ  + yQ) 


P + Qp  + JL  S1  + Js  S3  + R3 


(76) 


where  the  symbols  are  already  defined  for  Eqs.  (73)  - (75)- 

It  can  be  shown  that  the  major  component  of  the  imaginary  part  of  R^  arises 


from  J. 


V 


Both  this  coefficient  and  Jg  were  omitted  from  the  L-L  swerve* 


The  reduction  will • probably  be  expanded  to  handle  slowly  varying  spin  by  the 

V 2 

insertion  of  P f°r  v in  S^.  This  will  be  useful  for  finned  missiles 

over  the  portion  of  trajectory  for  which  v is  within  %%  of  steady  state This 
now  completes  our  discussion  of  the  individual  reduction  processes. 

• It  is  interesting  to  consider  the  relative  time  requirements  of  the  com- 
puting machines  for  the  various  problems.  This  comparison  is  provided  by  .Table  I, 
Although  the  Bell  Computer  is  rather  slow,  it  has  the  important  feature  of  being 
practically  error-free . 


1 The  advantages  of  numerical  integration  for  the  calculation  of  are  being 

considered.  This  is  feasible  when  the' problem  is  done  on  the  high  speed 
Ordvac . 
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TABLE  I 


Bell  Computer 

Time  in  Minutes 

1.  Drag  Reduction 

U unknowns 

15-2$ 

2 . Yaw  Drag  Reduction 

$ unknowns 

90-l$0 

3.  Roll  Differential  Corr. 

li  unknowns 

60-90 

k . Epicyclic  law  Diff.  Corr. 

6 unknowns 

7n  + kO 

5.  Epicyclic  Yaw  Diff.  Corr. 

8 unknowns 

9n  + 50 

6.  Epicyclic  Yaw  Diff . Corr. 

10  unknowns 

a . First  iteration 

12n  + 6$ 

b . Later  iteration 

lOn  + 6$ 

7.  Tricyclic  Yaw  Diff.  Corr. 

8 unknowns 

lUn  + 5$ 

8.  Tricyclic  Yaw  Diff.  Corr. 

10  unknowns 

a.  First  iteration 

l6n  + $0 

b . Later  iterations 

lUn  + 50 

9.  Epicyclic  Swerve 

6 unknowns 

6n  + 30 

10.  Epicyclic  Swerve  ' 

8 unknowns 

a.  Short 

7n  + 70 

b . Long 

12n  + 70 

11 . LF  Swerve 

6 unknowns 

a.  Short 

20n  + 30 

b.  Long 

2$n  + 30 

12 . LFS  Swerve 

7 unknowns 

a . Short 

23n  + UO 

b . Long 

28n  + UO 

13 . Tricyclic  Swerve 

10  unknowns 

6n  + UO 

lU,  Tricyclic  Swerve  (L-L) 

6 unknowns 

lOn  + 30 

* All  times  are  for  one  run  only,  n is  the  number  of  stations  in  the  round. 

**  Short  swerve  reductions  can  be  performed  only  when  the  yaw  reduction  has 
been  done  on  the  Bell  Computer, 
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TABLE  I (Con't) 


B Edvac 


Minutes 


1.  Epicyclic  Yaw 

2.  Epicyclic  Swerve 


10  unknowns 
8 unknowns 


2 •$ 

6.6 


Due  to  delays  of  tape  handling  the  usual  time  for  two  iterations  of  yaw 
plus  swerve  is  about  - 1 $ min. 


C Ordvac 


1.  Center  of  mass  reduction 

2 . Drag  Reduction 

3.  Roll  Diff . Com. 

U . Epicyclic  Yaw  Com. 
Epicyclic  Swerve  (L-F) 


U unknowns 
k unknowns 
10  unknowns 
6 unknowns 


(Un  + 30)  sec. 

2 min.  (max.) 

3 min. 

5 min. 

U min. 


t 

1 

I 


The  usual  time  for  two  iterations  of  yaw  plus  swerve  is  about  twelve  minutes. 


All  reductions  on  the  Ordvac  employ  the  binary  point  routines  described 
in  [12]  . 


}! 
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CRITERIA  FCR  QUALITY  OF  RESULTS 

In  conclusion  we  will  consider  the  various  ways  in  which  the  quality  of 
the  results  of  the  reductions  described  above  may  be  estimated.  There  are 
essentially  three  ways  by  means  of  which  this  can  be  done. 

First  the  sizes  of  the  residuals  of  each  fitted  curve  may  be  compared 
with  the  accuracy  of  the  corresponding  instrumentation.  This  comparison  will 
reveal  computational  errors  and  identify  rounds  for  which  the  particular  theory 
is  inappropriate  but  will  tell  very  little  about  the  accuracy  of  determination 
for  individual  coefficients  . 


The  second  method  describes  the 'internal  consistency  of  the  measurements 
and  involves  the.  computation  of  the  statistically  defined  standard  errors-*-  of 
each  aerodynamic  coefficient . Briefly,  if  X^  is  a parameter  of  one  of  the  least 


squares  fits  described  above,  it  can  be  shown  that  «Y,  the  standard  error  in  X,  , 

Ai  1 


may  be  computed  by  the  relation: 


where 


(standard  error  of  fit) 


(77) 


R residuals  of  fit 
n number  of  measurements 
r number  of  unknowns 


Ay,  elements  of  the  inverse  of  the  matrix  formed  by  the  coefficients  of 
the  least  squares  normal  equations . 

If  an  aerodynamic  coefficient  is  a function  f (X^,  X„)  of  two  of  these  parameters 
the  following  relation  is  used  to  calculate  the  coefficient's  standard  error: 

k\«  (70) 

£X„/  X,  X, 


-i  /di  _ \C  . /Oi  „ to i 

IVAi»  V ax1  xx  ax2  x2  \sxi: 


xi  2 


where 


*1*2 


A12e2'A21£2 


The  generalizations  of  Eq.  (78)  are  clear. 


1 The  commonly  used  probable  errors  may  be  calculated  by  multiplication  of 
the  standard  errors  by  ,67U5> . 


38 


I 


I 

j; 

» 

A 


I; 

V 


Experience  has  shown  that  the  relative  magnitudes  of  the  errors  for  the 
rounds  of  a program  are  quite  useful.  For  example  the  usual  size  of  errors 
in  Kjj  for  bodies  of  revolution  is  about  ,03.  Hounds  which  have  errors  in 

which  are  double  this  are  inspected  with  some  care.  The  absolute  magnitude 
of  statistical  errors  is  usually  much  smaller  than  the  round  to  round  scatter 
and  so  is  given  little  weight. 

The  third  and  most  effective  approach  examines  the  effect  which  the 
aerodynamic  coefficient  in  question  has  on  the  motion  and  determines  criteria 
for  the  motion  which  must  be  satisfied  in  order  that  this  effect  be  measurable . 
For  example,  in  order  to  determinate  the  drag  coefficient  the  observed  de- 
celeration of  the  missile  must  be  measurable  in  terms  of  timing  accuracies  of 
5/8  [i  sec and  distance  accuracy  of  either  ,0l5"  for  the  aerodynamics  range 
or  ,01*  for  the  transonic  range.  This  distinction  is  quite  important.  If  an 
extremely  heavy  missile  were  fired  with  no  yawing  motion  and  a drag  reduction 
performed,  the  residuals  of  fit  would  quite  probably  be  very  good.  The  stan- 
dard error  for  could  be  quite  reasonable  and  yet  the  would  be  meaning- 
less due  to  the  fact  that  the  actual  deceleration  would  be  smaller  than  could 
be  measured  accurately.  A second  example  could  be  constructed  from  the  roll 
reduction  of  a missile  whose  rolling  velocity  is  within  5$  of  steady  state 
throughout  its  observed  flight.  The  fit  could  be  excellent  and  the  calculated 
value  of  completely  worthless} 

Since  the  analyses  of  the  drag  and  roll  reductions  have  already  been  made 
in  [Ul  and  [lO]  respectively,  we  will  only  discuss  criteria  for  the  yaw  and 
swerve  reductions.  Experience  with  the  graphical  portion  of  the  yaw  reductions 
shows  us  the  and  K2Q  must  be  of  reasonable  size.  Since  depends  primarily 

on  ^ and  02  and  these  are  quite  well  determined  for  K^'s  over  five  times  the 

experimental  error  of  .001  rad,  the  criterion  for  K^.  is  that  K-^  5-  .005  and 

KgQ  ^ .005.  If  a spin  reduction  has  been  performed,  only  one  of  the  's  need 

be  larger  than  .005  while  0^-for  the  other  is  computed  from  Eq.  (36).  Since  II-  and  T 
depend  on  the  shrinking  or  grot-ring  of  the  two  arms,  the  requirement  on  the 
amplitudes  of  the  two  modes  must  be  strengthened.  Experience  has ' determined 
that  the  lower  bound  of  .007  for  the  amplitudes  of  both  modes  is  sufficient 
for  reasonable  determination  of  H and  T , 


The  criteria  for  swerve  reductions  requires  a little  more  algebraic  work. 
First  we  must  define  measures  of- the  effect  of  K^,  and  Kg  on  the  swerving 

motion.  In  Eq.  (67)  only  Rg  need  be  considered  for  spinning  bodies  of 

revolution.  For  this  case  the  lift  swerve  (R^^  anc*  Magnus  swerve  (R2 )p 

are  defined  by  the  relations: 


“ JL  K20  y*2I  + a222 
(R2)f-JFK20  /b212  + b222 


(79) 

(80) 


f 
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We  then  require  tl  at^)^  or  (Rg)^  be  five  times  the  experimental  accuracy  in 

calibers  in  order  that  the  corresponding  coefficient  be  determined . (This 
transverse  distance  accuracy  for  the  Transonic  Range  is  .01  ft,  while  for  the 
Aerodynamics  Range  it  is  .015  in.)  For  the  six  unknowns  swerve  described  by 
Eqs . (6?)  and  (70),  the  following  definitions  are  natural: 


JL  Jho 


h K Js 


' Z 2 ' 7 Z 2 

ho  ho 


2 2 

°21  + c22 


Even  simpler  definitions  are  possible  for  Eq.  (73) 


SL  ‘ JL 


( sx|  2 ) 


average 


asi|  2 


SF  ‘ vo  JF 


Ss-Js 


'(|S,  t 


average 


( s,  2) 


average 


= y0  jf 

' n 

* JsVwr 


where  summations  are  over  all  stations  and  n is  the  number  of  stations'1".  The 
criteria  for  (R)^,  (R)s,  Sp,  Ss  is  the  same  as  that  for  (Rg)^  (R2)p. 

The  tricyclic  swerve  may  be  treated  in  the  same  manner. 


1 Since  2 appear  as  a coefficient  in  the  normal  equations,  definitions 

(83)  - (85)  are  especially  convenient. 
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TABLE  OF  SYMBOLS 

Axial  moment  of  inertia $ constant  in  Roll  Eq.  (31) 

Element  of  inverse  of  normal  matrix  for  least  squares  [’Eq.  (77)] 
Transverse  moment  of  inertiaj  constant  in  Roll  [Eq.  (31)] 
Constant  in  modified  roll  equation  [Eq.  (31*)] 


D * JD  ‘ 4-  JA 
A 

D * D at  M * M. 


E Constant  in  modified  roll  equation  [Eq.  (31’)] 

(F^^jF^)  Aerodynamic  Force  vector  in  missile  coordinate  system 

(fx,fy>f2)  Aerodynamic  Force  vector  in  range  coordinate  system 

a . iS|  [j „ ♦ kj2  JH  - Iv] 

U 2 

H * JL  ‘ JD  * IT  JH 

I1(p)  Integral  defined  for  Eq.  (26) 

Ip(p)  Integral  defined  for  Eq.  (26) 


■ K-§>  \ - *2'2  v]  \ 


Spin  deceleration  moment  coefficient 
Spin  producing  moment  coefficient 
Drag  force  coefficient  (trajectory  drag) 
Axial  drag  force  coefficient 


Yaw  drag  coefficient 

6 a62 


Slope  of  Mach  number  drag  curve 


Magnus  force  coefficient 


1 A few  symbols  unfortunately  have  two  meanings.  These  are  separated  by  semi- 
colons. The  definition  which  applies  may  easily  be  determined  by  the  context* 


bZ 


n * 


kh 

*l 

*m 

F-N 

*s 

Kp 

*XF 

^XT 

Kio 


TABLE  OF  SYMBOLS  (Con't) 
Damping  moment  coefficient 
Lift  force  coefficient 
Overturning  (righting)  moment  coefficient 
Asymmetrical  moment  coefficient 

Normal  force  coefficient 
Asymmetrical  normal  force  coefficient 

Damping  force  coefficient 

Magnus  moment  coefficient 

Magnus  force  coefficient  due  to  cross  spin 

Magnus  moment  coefficient  due  to  cross  spin 

Amplitude  of  the  ith  mode  of  yawing  motion 


i 0. 


K.  =*  K.  e 'xo 
i xo 


Distance  from  missile* s reference  point  to  center  of  mass 


-2. 


rW 


M = kg  J^j  or  Mach  number 
(M^MgjM^)  Aerodynamic  moment  vector 
P Constant  in  Eq.  (f?8)j  Pressure 

Vapor  pressure  of  water 
Pressure  of  dry  air 
Constant  in  Eq,  (£8) 

Residuals 

Expression  defined  for  Eq.  (67) 

Ri  * Ril  + **12 

(Rg)^,(R)^  Radii  of  lift  swerve  motion  for  Eqs . (67)  and  (69)  respectively 
(Rg)p  Radius  of  Magnus  swerving  motion  for  Eq,  (67) 


R 


R.  . 
10 
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TABLE  OF  SYMBOLS  (Con't) 


(E)s 
®X>  % 


Radius  of  damping  force  svr:*-nring  motion  for  Eq.  (6?) 
Constants  in  Eqs  . (6)  and  (7) 

Quantities  appearing  in  Eq.  (73) 

Radius  of  Magnus  swerving  motion  [Eq.  (8l;)3 

Radius  of  lift  swerving  motion  [Eq.(83)3 

Radius  of  damping  force  swerving  motion  [Eq . (85)3 


_2 

T = - kj_  J,pj  Absolute  temperature 

Vs  Velocity  of  sound 

VSD  Velocity  of  sound  in  dry  air 


a,b 


iJ 

b. , 
ij 

(V°2> 

(OpCj) 

c.  , 

10 


Constants  defined  by  D*(a  + bM)  ^ 
Coriolis  acceleration 
Coefficient  of  drag  equation  [Eq.  (21)3 
Quantities  defined  for  Eq.  (67) 
Quantities  defined  for  Eq..  (67) 
x and  y coordinates  of  horizontal  spark 
x and  y coordinates  of  vertical  spark 
Quantities  defined  for  Eq.  (67) 
Diameter 

Quantities  defined  for  Eq.  (67) 


Unit  vectors  along  1,2,3  axis  on  earth's  surface 


Unit  vectors  along  x,y, z axi3  located  on  the  range 


g 

k 


Acceleration  due  to  gravity 
Axial  radius  of  gyration  in  calibers 


hh 


TABLE  OF  SYMBOIS  (Con't) 


k9  3 /~*  Transverse  radius  of  gyration  in  calibers 

Vr.d 

m Lass 

in  = v2  - UM  - H2 

(m,n,p)  Direction  cosines  of  missile's  axis 

n Number  of  observations 


P’d 


u 

U1 

x,y,z 


Number  of  unknowns;  quantity  defined  in  Eq.  (29) 
Steady  state  spin 

_p 

Constants  defined  by  k-^  * - D(sq  + 3*P) 

0 

Missile's  velocity 

Axial  component  of  missile's  velocity 
Positional  coordinates  in  range  system 

(XjpO,Zjj),  (0,yy,Zy)  Coordinates  of  plate  measurements 

(WR)Zn)  Space  coordinate  of  reference  point  on  missile 
^xcm,ycm,zcm^  sPace  coordinate  of  center  of  mass  of  missile 
x^y^  x and  y components  of  drift 

x^,y^  x and  y components  of  Coriolis  deflection 

yQ  Gravity  drop 


on  Exponential  damping  coefficient  of  ith  mode 

p Exponent  defined  for  Eq.  (71) 

2 

h m - v + i(^op+  K f) 

A A 

p±  - P±  + Dp 

-Jn 

T “ I+a7bM 

o 

r - 02o  " ^io  ’27  ^Eq* 
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TABLE  OF  SYMBOLS  (Can't) 


6 = | X 1 magnitude  of  yaw 
6^  Mean  squared  magr 


n - r, 


Mean  squared  magnitude  of  yaw 
Standard  error  of  fit 


Standard  error  in  i 


Perturbation  term  arising  from  slowly  varying  spin  fEq.  (35)] 
~>1  - arccos « - (ql  * fl2)  [Eq.  (29)] 


0 Roll  angle;  latitude 

X « Xjj  + i Xy  = Xg  + i X^ 

Xr,  Xy,  Xg,  X^  Components  of  yaw  in  radians  along  horizontal,  vertical 
J 1 and  2 axis  respectively 

Xe  Magnitude  of  asymmetry  angle 


Xr  Yaw  of  repose 

» - “2d  + 1 “j*  nondi 

v * ^1^  nondi 


nondimens ional  cross  angular  velocity 


nondimensional  spin 


- A 

V “S  v 
P * Pa  + Pw 


air  density 


Pa,  pw  "partial  densities"  defined  after  Eq.  (2) 

$i  * ^ jy-  orientation  of  i,  the  mode,  in  complex  yaw  plane 


Azimuth 


Magnitude  of  earth* s angular  velocity 


Missile's  angular  velocity 


( ) Derivative  with  respect  to  time 

t 

( ) ’ Derivative  with  respect  to  p 

\ 7 Average  or  conjugate  of  complex  quantity 


1*6 
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